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Abstract. Consider a random multigraph G* with given vertex de- 
grees di, . . . , d„, contracted by the configuration model. We show that, 
asymptotically for a sequence of such multigraphs with the number of 
edges I — > oo, the probability that the multigraph is simple stays 
away from if and only if J2^d] = 0{J2^di). This was previously 
known only under extra assumtions on the maximum degree maxi di . 
We also give an asymptotic formula for this probability, extending pre- 
vious results by several authors. 



1. Introduction 

If n > 1 and (dj)" is a sequence of non-negative integers, we let G{n, (dj)") 
be the random (simple) graph with the n vertices 1, . . . , n, and with vertex 
degrees di, . . . , d„, uniformly chosen among all such graphs (provided that 
there are any such graphs at all; in particular, t° be even). A 

standard method to study G{n, {di)^) is to consider the related random 
multigraph G* (n, (di)" ) defined by taking a set of di half-edges at each vertex 
i and then joining the half-edges into edges by taking a random partition of 
the set of all half-edges into pairs; see Section 2 for details. This is known as 
the configuration model, and such a partition of the half-edges is known as a 
configuration; this was introduced by Bollobas [2], see also Section II. 4 of [3]. 
(See Bender and Canfield [1] and Wormald [14, 15] for related arguments.) 

Note that G*{n, (di)") is defined for all n > 1 and all sequences (dj)" such 
that dj is even (we tacitly assume this throughout the paper), and that 
we obtain G{n, (dj)") if we condition G*{n, (d^)^) on being a simple graph. 
The idea of using the configuration method to study G(n, (dj)") is that 
G*{n, {di)i) in many respects is a simpler object than G{n, (d^)"); thus it is 
often possible to show results for G{n, {di)^) by first studying G*{n, {di)^) 
and then conditioning on this multigraph being simple. It is then of crucial 
importance to be able to estimate the probability that G*{n, (d^)") is simple, 
and in particular to decide whether 

liminf P(G*(n, (di)^) is simple) > (1.1) 
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for given sequences {di)i = (c?!"^)" (depending on n > 1). (Note that (1.1) 
implies that any statement holding for G*{n, (d^)") with probability tending 
to 1 does so for G{n, ((i;)i) too.) 

A natural condition that has been used by several authors using the con- 
figuration method (including myself [7]) as a sufficient condition for (1.1) 
is 

n n 

Y^di = Q{n) and ^df = Oin) (1.2) 
i=l i=l 

together with some bound on maxjdj. (Recall that A = Q{B) means that 
both A = 0(B) and B = 0(A) hold.) Results showing, or implying, that 
(1.2) and a condition on maxj di imply (1.1) have also been given by several 
authors, for example Bender and Canfield [1] with maxjdj = 0(1); Bollobas 
[2], see also Section II. 4 in [3], with maxj < \/2 log n — 1; McKay [10] 
with maxjdj = o(n^/^); McKay and Wormald [13] with maxjdj = o(n^/^). 
(Similar results have also been proved for bipartite graphs [9], digraphs [5], 
and hypergraphs [4].) 

Indeed, it is not difficult to sec that the method used by Bollobas [2, 3] 
works, assuming (1.2), provided only maxjdj = o(n^/^), see Section 7. This 
has undoubtedly been noted by several experts, but we have not been able 
to find a reference to it in print when we have needed one. 

One of our main result is that, in fact, (1.2) is sufficient for (1.1) without 
any assumption on maxj di, even in cases where the Poisson approximation 
fails. Moreover, (1.2) is essentially necessary. 

We remark that several papers (including several of the references given 
above) study P(G*(n, (d^)") is simple) from another point of view, namely 
by studying the number of simple graphs with given degree sequence (c/^)". 
It is easy to count configurations, and it follows that this number equals, 
with A'^ the number of edges, see (1.3) below, 

^^^^^P(G*(n, (d,)?) is simple); 

such results are thus equivalent to results for P(G* (ra, (c?i)i ) is simple) . How- 
ever, in this setting it is also interesting to obtain detailed asymptotics when 
P(G*(n, {di)i) is simple) — > 0; such results are included in several of the ref- 
erences above, but will not be treated here. 

We will throughout the paper let be the number of edges in G*{n, {di)"^). 
Thus 

n 

2N = J2di- (1-3) 

i=l 

It turns out that it is more natural to state our results in terms of N than 
n (the number of vertices). We can state our first result as follows; we use 
an index v to emphasize that the result is asymptotic, and thus should be 
stated for a sequence (or another family) of multigraphs. 
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Theorem 1.1. Consider a sequence of random multigraphs G* = G* {uy, {d['^^)i) . 

Let Ny = i Ylii d'f\ tf^^ number of edges in G*, and assume that, as u ^ oo, 
N„^oo. Then 

(i) liminf^^ooF(G* is simple) > if and only ifY^Mf^f = 0{Ny); 
in) \im.y^^T{G* is simple) = if and only if^i{d['^^)'^/Ni, — ^ oo. 

In the sequel we will for simplicity omit the index v, but all results should 
be interpreted in the same way as Theorem 1.1. 

Usually, one studies G*{n,{di)i) as indexed by n. We then have the 
following special case of Theorem 1.1, which includes the claim above that 
(1.2) is sufficient for (1.1). 

Corollary 1.2. Let {di)1 = (df ^)? be given for n > 1. Assume that N = 
Q{n). Then, as n ^ oo, 

(i) liminf„_>ooP(G*(n, (c/j)") is simple) > if and only if ^i{d\'^^)'^ = 
0{n), 

(ii) P(G*(ra, {di)i) is simple) if and only if ^iidf^)"^ /n —>■ oo. 

Remark 1.3. Although we have stated Corollary 1.2 as a special case of 
Theorem 1.1 with N = 0{n), it is essentially equivalent to Theorem 1.1. In 
fact, we may ignore all vertices of degree 0; thus we may assume that di > 1 
for all i, and hence 2N > n. If further J2i di = 0{N), the Cauchy-Schwarz 
inequality yields 

2N = ^di<(n^df^ =0{V^), 
i=i i=i 

and thus = B(n). In the case J2i df/n^ oo, it is possible to reduce some 
di to 1 such that then N = \ Y^idi = 6(n) and still J2id'i/N — oo; we 
omit the details since our proof does not use this route. 

Our second main result is an asymptotic formula for the probability that 
G*{n, idi)i) is simple. 

Theorem 1.4. Consider G* in, [d-i)'^) and assume that N :=^J2i di oo. 
Let \ij := -yj diidi — l)dj{dj — 1)/(2N); in particular Xa = di{di — 1)/{2N). 
Then 

F{G*{n,{di)'l) is simple) 

= exp(-i^Aii-^(A,j-log(l + A,j))) +o(l); (1.4) 

i i<j 
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equivalently, 

F{G*{n,{dJ{) is simple) 

+ o(l). (1.5) 

In many cases, ^jC?f(c?i — 1)^ in (1.5) may be replaced by the simpler 
'^j^dj; for example, this can be done whenever (1.1) holds, by Theorem 1.1 
and (2.4). Note, however, that this is not always possible; a trivial counter 
example is obtained with n = 1 and di = 2N oo. 

In the case maxj di = o{N^/'^), Theorem 1.4 simplifies as follows; see also 
Section 7. 

Corollary 1.5. Assume that N ^ oo and max,; d, = o{N^/^). Let 

A:= (1.6) 

i=l 

Then 



P(G*(n, (di)^) is simple) = exp(-A - A^) + o(l) 



A\ 2N J 4. 

This formula is well known, at least under stronger conditions on maxj dj, 
see, for example. Bender and Canfield [1], Bollobas [3, Theorem 11.16], 
McKay [10] and McKay and Wormald [13, Lemma 5.1). 

2. Preliminaries 
We introduce some more notation. 

We will often write G* for the random multigraph G*{n, (dj)"). 

Let Vn = {1, ■ ■ ■ ,n}; this is the vertex set of G*{n, (dj)"). We will in the 
sequel denote elements of Vn hy u,v,w, possibly with indices. Vn is also the 
vertex set of the complete graph K^, and we let En denote the edge set of 
Kn] thus En consists of the (2) unordered pairs {v,w}, with v,w EVn and 
V ^ w. We will use the notation vw for the edge {v,w} G En- 

For any multigraph G with vertex set Vn, and u G Vn, we let Xu{G) be 
the number of loops at u. Similarly, if e = vw G En, we let X(.(G) = X^^(G) 
be the number of edges between v and w. We define further the indicators 

Iu{G) := l[Xu{G) > 1], ueVn, 

JeiG):=l[XeiG)>2], eeEn, 

and their sum 

Y{G) := ^ Iu{G) + ^ MG). (2.1) 
uev„ eeEn 
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Thus G is a simple graph if and only iiY(G) = 0, and our task is to estimate 
P(y(G*) = 0). 

As said above, the idea of the configuration model is that we fix a set of 
dy half-edges for every vertex v; we denote these half-edges by v^^\ . . . , v^'^^\ 
and say that they belong to or are at v. These sets are assumed to be 
disjoint, so the total number of half-edges is dy = 2N. A configuration is 
a partition of the 2N half-edges into N pairs, and each configuration defines 
a multigraph with vertex set Vn and vertex degrees dy by letting every pair 
{x, y} of half-edges in the configuration define an edge; if a; is a half-edge at 
V and y is a half-edge at w, we form an edge between v and w (and thus a 
loop ii V = w). We express this construction by saying that we join the two 
half-edges x and y to an edge; we may denote this edge by xy. Recall that 
G* is the random multigraph obtain from a (uniform) random configuration 
by this construction. 

We will until Section 6 assume that 

Y,dl = 0{N), (2.2) 

V 

i.e., that J2ydl < CN for some constant C. (The constants implicit in the 
estimates below may depend on this constant C.) Note that an immediate 
consequence is 

max dy = 0{N^/^) = o{N). (2.3) 

V 

We may thus assume that N is so large that vnscy^^dy < A^/10, say, and 
thus all terms like N — dy are of order A^. (The estimates we will prove are 
trivially true for any finite number of N by taking the implicit constants 
large enough; thus it suffices to prove them for large N.) 

Note further that (2.2) implies, using (2.3), that for any fixed k >2 

< (ma^d,)'=-2^d2 = 0{N''/^). (2.4) 

V V 

We further note that we can assume > 1 for all v, since vertices with 
degree may be removed without any difference to our results. (This is 
really not necessary, but it means that we do not even have to think about, 
for example, d~^ in some formulas below.) 

We will repeatedly use the subsubsequence principle, which says that if 
{xn)n is a sequence of real numbers and a is a number such that every 
subsequence of {xn)n has a subsequence that converges to a, then the full 
sequence converges to a. (This holds in any topological space.) 

We denote the falling factorials by x- := x{x — 1) ■ ■ ■ {x — k + 1). 

3. Two PROBABILISTIC LEMMAS 

We will use two simple probabilistic lemmas. The first is (at least part 
(i)) a standard extension of the inclusion-exclusion principle; we include a 
proof for completeness. 
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Lemma 3.1. Let W he a non-negative integer-valued random variable such 
thatER^ < oo for some R>2. 

(i) Then, for every j > 0, 

k=j 

(ii) More generally, for every random variable Z such that M{\Z\R^) < 
oo for some R> 2, and every j > 0, 

K{Z . 1[W = j]) = f^i-l)'-^ Q ^ E{ZW^). 

Proof For (i), let f{t) := E{t^) = = i)P be the probability 

generating function of W\ this is by assumption convergent for |t| < R, at 
least. If |t| < R—\ we have 

/(i+i) = E(i+i)^ = f;Ef^)t^ 



and thus, if \t\ < R — 2, 



/(*) = E^U = E^(^-AOE •r'(-i)'"'- (3.1) 

t-n V / i.-n ,-n \J / 



k=0 ^ ' k=0 j=0 

The double series is absolutely convergent since 

f; f; E(l^Vfc!) |ip = E E {\t\ + 1)^= = f{\t\ + 2) < oo. 

fc=Oi=0 k=0 ^'^^ 

Hence the result follows by extracting the coefficients of P in (3.1) 

Part (ii) is proved it the same way, using instead f{t) := E{Zt^). □ 

The next lemma could be proved by Lemma 3.1 if made the hypothesis 
somewhat stronger, but we prefer another proof. 

Lemma 3.2. Let {Wy)i, and {W^)i, be two sequences of non-negative integer- 
valued random variables such that, for some R> 1 

supE(i2^^'''^) < oo (3.2) 

and, for each fixed k > 1, 

E{W^) - E{W^) ^0 as u^oo. (3.3) 

Then, as u ^ oo, 

F{W^ = 0) - F{W„ = 0) ^ 0. (3.4) 
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Proof. By the subsubsequence principle, it suffices to prove that every sub- 
sequence has a subsequence along which (3.4) holds. Since (3.2) implies 
that the sequence {Wu)u is tight, we can by selecting a suitable subsequence 

assume that Wy — ^ W for some random variable W (see Sections 5.8.2 and 
5.8.3 in Gut [6]). Moreover, (3.2) implies uniform integrability of the powers 
for each k, and we thus have, as — oo along the selected subsequence, 
E(W^^) ^ W.{W^) for every k and thus also E{W^) E(W^) (see Theo- 
rems 5.4.2 and 5.5.9 in [6]). By (3.3), this yields also E(TyJ^) ^ E(W^). 
Furthermore, (3.2) implies by Fatou's lemma (Theorem 5.5.8 in [6]) that 
E(i?'^) < liminfE(i2^'-) < oo, or E(e*^) < oo with t = \ogR > 0; hence 
the distribution of W is determined by its moments (see Section 4.10 in [6]). 
Consequently, by the method of moments (Theorem 6.7 in [8]), still along 

the subsequence, W and thus 

P(VF^ = 0) - P(W^^ = 0) ^ ¥{W = 0) - ¥{W = 0) = 0. □ 

Remark 3.3. The same proof gives the stronger statement 

dTv(W„ Wl) := I = j) - ^i^- = i)l ^ 0- 

i 

4. Individual probabilities 

We begin by estimating the probabilities P(/„(G*) = 1) and P(J^„,(G*) = 
1). The following form will be convenient. 

Lemma 4.1. Suppose ^^dl = 0{N). Then, for G* , and for all u,v,w G 
Vn, if N is so large that du < N, 

- logm = 0) = - logP(X„ = 0) = + O (^) 

and, with X^w := ^/dy{dv — l)dw{dw — 1)/{2N) as in Theorem I.4, 
-logP(J^^ = 0) = -logP(X^^ < 1) 

= - l0g(l + A„^) +Xyw+Ol ^ ^3 1 ■ 

Proof. The calculation for loops is simple. We construct the random con- 
figuration by first choosing partners to the half-edges at u, one by one. A 
simple counting shows that 

p(x.=o)=n(i-^^Az^) 

and thus, for large N, using — log(l — x) = x + 0{x^) when \x\ < 1/2, 

— (x„^o)^|(^.o(^)).*^.o(<). 
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For multiple edges, a similar direct approach would be much more com- 
plicated because of the possibility of loops aX v or w. We instead use 
Lemma 3.1(i), with W = Xy^- Wc may assumG du^di^j ^ 2, since tli6 

k 

result otherwise is trivial. X^w is the number of ordered fe-tuples of edges 
between v and w; the corresponding pairs in the configuration may be cho- 
sen in dydnj ways, and each such set of k pairs appears in the configuration 
with probability {{2N - 1){2N - 3) ■ ■ ■ {2N - 2k + 1))-^ Thus 



Ijj, dydy] dudyj T— T {dy l) {dyj l) 

~ {2N -!)■■■ (27V -2k + l) ~ 2^{N-l/2)h. ~ ii 2{N-l/2-i) ' 

(4.1) 

In particular, 

EXf^ = AL(l + 0(l/Ar)) (4.2) 
and if A; > 3, uniformly in k, 

fc-i 

^^^'2{N-l/2-i) 



= (1 + 0(l/Ar))AL n ^,t/v-!/9-,? - (4-^) 



Since d^ < N (for large N at least, see (2.3)), the ratios [dy — i)/{N 
1/2 — i) decrease as i increases; hence, for large and i > 2, 



dv-i ^ lit, - 2 ^ dy-l ^ \/dy{dy - 1) 



2{N -i-l/2) - 2N 2N 2N 
and (4.3) yields, uniformly for A; > 3, 

EXi < (1 + 0(l/iV))AL n ^--^ = (1 + 0(VA^))AL- (4.4) 
In the opposite direction, still uniformly for > 3, 

"^^Ww/Kw > ii > ii j2 

1=2 1=2 
fe— 1 fc— 1 ,2 1,2 

= J](l - iK)(l - i/d^) > 1 - Y,{i/dy + > 1 - - 

i=2 i=2 

Together with (4.4), this shows that, uniformly for > 3, 

EX,\, = A^^ (1 + O (A;2(d;i + d-i))) . (4.5) 

We now use Lemma 3.1(i), noting that trivially ¥,R^^^ < R'^^ < oo for 
every R. Thus, using (4.2) and (4.5) and observing that A^^ = 0(1) by 
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(2.3), 

P(X^^ = 0) = l-EX^^ + iAL(l + 0{N-^)) 

k=3 

= 1 - EX„^ + e-^"- - 1 + A^^ + O (A^^iV-i + Xl^e^^-{d-^ + rf-^)) 
= e-^"- + A„^ - EX^^ + O ((d^^ + dldl)N-^) . 
Similarly, still by Lemma 3.1(i), 
P(X„^ = 1) = EX„^ - AL(1 + OiN-')) 

= EX„,^. + A^^(e-^"» - 1) + O (aLa^^^ + Xiy^^id-^ + d"^)'^ 

= EX„^ + A^^e"^"™ - A^^ + O ((^^4 + dldl)N-^) . 

Summing these two equations we find 

P(X„, <!) = (! + X,^)e-^^- + O {{dldl + dldl)N-^) 

and the result follows, noting that (1 + Xyyj)e~^'"" is bounded below since 
A^^ = 0(1). □ 

5. Joint probabilities 

Our goal is to show that the indicators and Je{G*) are almost 

independent for different u and e; this is made precise in the following lemma. 
We define for convenience, for u G 14 and e = vw e En, 

fJ-u ■= dl/N and fXe ■= d„dw/N. (5.1) 

It follows easily from (4.1) and a similar calculation for loops that 

E(X„(G*)^) < ^li and E(Xe(G*)^) < i4, k> 1. (5.2) 
In particular, omitting the argument G*, 

P(/„ = l)=E/„<EX„</x„, 

F{Je = 1) = E Je < EX| < /X^. 



(5.3) 



More precisely, it follows easily from Lemma 4.1 that (for N large at least) 
P(4 = 1) = e(/Xu) and P(Je = 1) = 0(/«e) provided dy,dw > 2; this may 
help understanding our estimates but will not be used below. 

Lemma 5.1. Suppose ^^d^ = 0{N). Let I > and m > be fixed. For 
any sequences of distinct vertices ui, . . . ,ui E Vn and edges ei, . . . , G E^, 
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let ej = VjWj and let F be the set of vertices that appear at least twice in the 
list ui,... ,ui,vi,wi, . . .,Vm,Wm- Then, 

(I m \ I m 

n (G*) n Je, (G*) = n (g*)) n (g*)) 
1=1 j=i ^ i=i j=i 

(I m \ 

(N-^ + Y.dv')\{l^u.\{l^l). (5.4) 
veF 1=1 j=i ^ 

The implicit constant in the error term may depend on I and m but not on 
{ui)i and {ej)j. Ah similar statements below are to be interpreted similarly. 

The proof of Lemma 5.1 is long, and contains several other lemmas. The 
idea of the proof is to use induction in Z + m. In the inductive step we 
select one of the indicators, Je^ say, and then show that the product of the 
other indicators is almost independent of Xg^ , and thus of J^^ . In order to 
do so, we would like to condition on the value of Xg^. But the effects of 
conditioning on Xf,^ = k are complicated and we find it difficult to argue 
directly with these conditionings (see Remark 5.7). Therefore, we begin 
with another, related but technically much simpler conditioning. 

Fix two distinct vertices v and w. For < A; < min(d^,(i^), let Sk be 
the event that the random configuration contains the k pairs of half-edges 
{v^^\w^^^}, 1 < i < fc, and let the corresponding random multigraph, i.e., 
G* conditioned on S^, be denoted G^. G*^ thus contains at least k edges 
between v and w, but there may be more. Note that Gq = G*. 

We begin with an estimate related to Lemma 5.1, but cruder. 

Lemma 5.2. Suppose = 0{N). Let I, m and ri, . . . ,ri, si, . . . , s^a be 

fixed non-negative integers. For any sequences of distinct vertices ui, . . . ,ui G 
Vn and edges ei, . . . , G En, 

y I m \ ^ I m \ 

H=l j=l ^ \=\ j=l ' 

In particular, 

(I rn \ /I m \ 

n4,(G*)n<^e,(G*) =o n^n^n^ • (^-e) 
i=i j=i ^ \=i j=i 

The estimates (5.5) and (5.6) hold with G* replaced by too, uniformly 
in k, provided the edges ei,...,em are distinct from the edge vw used to 
define G\. If vw equals some Cj, then (5.5) still holds for G*j^, if we replace 

X(,j by Xg. — k when ej = vw. 

Proof. We argue as for (4.1). Let, again, ej = VjWj and let t := ri + • • • + 
n + si + ■ ■ ■ + Sm- The expectation in (5.5) is the number of t-tuples of 
disjoint pairs of half-edges such that the first ri pairs have both half-edges 
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belonging to ui, and so on, until the last Sm that each consist of one half- 
edge at Vm and one at Wm, times the probability that a given such t-tuple 
is contained in a random configuration. The number of such f-tuples is at 
most Yli=i W%i{dv,du,^Y' and the probability is ((2iV - 1) • • • {2N - 
2t + 1))~^ < AT"* (provided N > 2t). The estimate (5.5) follows, recalling 
(5.1), and (5.6) is an immediate consequence since 7„ < X„ and Je < Xe- 

The same argument proves the estimates for G^,. There is a minor change 
in the probability above, replacing N hy N — 2k] nevertheless, the estimates 
are uniform in k because k < dy = O^N^/"^). (There may also be some t- 
tuples that are excluded because they clash with the special pairs {v^'^^ , u;^') }, 
i = 1, . . . , A;; this only helps.) □ 

Let . . . ,ui G Vn and ei , . . . , e En be as in Lemma 5.1, and assume 
that m > 1. We choose v = vi and w = wi, so ei = vw, for the definition 
ofG*. 

If A; > 1, we can couple G^. and G^_^ as follows. Start with a random con- 
figuration containing the k special pairs {v^'^\w^^^}. Then select, at random, 
a half-edge x among all half-edges except v^^\ . . . , v^''\w^'^\ . . . , ■w^''~^\ If 
X = w'^^^ do nothing. Otherwise, let y be the half-edge paired to x\ remove 
the two pairs {v^^\w'^^^} and {x,y} from the configuration and replace 
them by {v^^\x^ and {w^^\y^. (This is called a switching; see McKay [10] 
and McKay and Wormald [11, 13] for different but related arguments with 
switchings.) 

It is clear that this gives a configuration in £k-i with the correct uniform 
distribution. Passing to the multigraphs, we thus obtain a coupling of G^. 
and G*i^_^ such that the two multigraphs differ (if at all) in that one edge 
between v and w and one other edge have been deleted, and two new edges 
are added, one at v and one at w. 

Let Z denote the product 171=1 ^^=2 "^ej of the chosen indicators ex- 
cept Je^. Define Fi C {v,w} to be the set of endpoints of vw = e\ that also 
appear as some Ui or as an end-point of some other e^; thus F\ = Fn {f , vj). 
We claim the following. 

Lemma 5.3. Suppose = 0{N). With notations as above, uniformly 

in k with 1 < k < mm{dy,dw), 

(I m \ 

(N-'+Y,d-')ll^u^l[filX (5.7) 
v£Fi i=l j=2 ' 

Proof. We use the coupling above. Recall that Z = or 1, so if Z{G'^) and 
Z{Gl_^) differ, then one of them equals and the other equals 1. 

First, if Z{Gl) = 1 and Z{Gl^^) = 0, then the edge xy deleted from Gl 
must be either the only loop at some Ui, or one of exactly two edges between 
Vj and Wj for some j > 2. Hence, for any configuration with Z{G*^,) = 1, 
there are less than I + 2m such edges, and the probability that one of them 
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is destroyed is less than {I + 2m) /{N - k) = 0{l/N). Hence, 

^{Z{Gl) > Z{Gl_,)) = 0{nZ{Gl))/N). (5.8) 

Define M := l\\=i M«i 11^2 /^e, - By Lemma 5.2, EZ(Gp = 0(M), so the 
probabihty in (5.8) is 0{M/N), which is dominated by the right-hand side 
of (5.7). 

In the opposite direction, Z{Gl.) = and Z{G^_^) = 1 may happen in 
several ways. We list the possibihties as fohows. (It is necessary but not 
necessarily sufficient for Z^G^.) < Z{Gl._i) that one of them holds.) 

(i) V is an endpoint of one of the edges 62, . . . , e^, say v = V2 so 62 = 
VW2; the new edge from v goes to W2; there already is (exactly) one 
edge between v and W2 in G^; if we write Z' = ni=i Hjla "^ej > so 
that Z = Je^Z', then Z'(G^) = 1. 

(ii) V equals one of ui,. . . ,ui, say ^ = 1x1; the new edge from i; is a 
loop; if we write Z' = 11^=2 -^wi 11^2 "^e^ ' so that Z = lu^Z', then 
Z'iGl) = 1. 

(iii) Two similar cases with v replaced by w. 

(iv) Both V and w are endpoints of edges ej, say v = V2 and w = W3, 
so that 62 = VW2 and 63 = WV3; the two new edges go from v 
to W2 and from w to V3; there are already such edges in G^; if 
Z" = n!=i4in7=4<^e,> so that Z = Je,JesZ", then Z"(Gl) = 1, 
where Gj. is Gj^ with one edge between W2 and vs deleted. 

(v) Both V and w equal some Ui, say v = ui and w = U2; the new 
edges from v and w are loops; if Z" = n!=3 -^wi 11^2 "^e^ > so that 
Z = luJu^Z", then Z"(G^) = 1. 

(vi) A similar mixed case where, say v = V2 and w = ui. 

(vii) The same with v and w interchanged. 

Consider case (i). For any configuration, the probability that the new 
edge from v goes to W2 is dw^/{2N — 2k + 1) = 0{dw2/N)- Since we also 
need Z'{G*j^) = 1 and Xe2(G^) > 1, the probability of case (i) is at most 

0{d^,N-^W.{X,,{Gl)Z'{Gl))). 

Now, by Lemma 5.2, for convenience omitting the arguments here and 
often below in this proof, 

(I ra \ / I m \ 

i=l j=3 ^ \=1 i=3 ^ 

Moreover, dy,^/N = He^/d^, so the probability of case (i) is 0{M/dy); note 
that the case only can happen if v G Fi , so this is covered by the right-hand 
side of (5.7). 

Case (ii) is similar (and slightly simpler). 

Case (iii) occurs, by symmetry, with probability 0(M/dw), and only if 
w e Fi. 
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In case (iv), the other destroyed edge must go between 102 and ^3. For any 
configuration, the probabihty that such an edge is chosen is 0{Xyj^v^/N). 
We study two subcases. If one of the edges 64,..., equals W2V^^ say 
64 = W2V'i, then we, moreover, need at least three edges between W2 and 
vz in Gj^, since one of them is destroyed. We also need Xe^{G\) > 1 and 
Xe^iG].) > 1. Thus the probability of this case then is 

0(iV-lE(X^,,3Xe,Xe3Z"(G;))) =0(iV-lE(Xe,Xe3Xe,l[Xe, >3]Z")). 

By Lemma 5.2 we have 

(I m \ 

l[x^^.x,,x,,xlllxf.] 
i=l j=5 ' 

(I m 
i=i 7=5 



j=5 

= 0{MHeJ{l^e2l^e3))- 

In this case we have /Ue2A*e3 = AteiA*e4; so the probability is 0{N~^M / jif,^) = 
0{M/d^d^). 

In the second subcase, W2V^ does not equal any of 64, ... , e^. We then 
obtain similarly the probability 



(I m 
i=l j=4 

= 0(iV"^MyU^2^3/(/Xe2Me3)), 

which again equals 0{M / {dydw)). Finally, note that in case (iv), Fi = 
{v,w}. 

In case (v), the other destroyed edge is also an edge between v and w] given 
a configuration, the probability of this is 0((X„^ — k)/N). The probability 
of case (v) is thus 

(I m 
i=3 j=l 

= 0{M/{dyd^)). 

Fi = {v, w} in case (v) too. 

Cases (vi) and (vii) are similar to case (iv), and lead to the same estimate. 
Again Fi = {v,w}. 

By (5.8) and our estimates for the different cases above, the probability 
that Z{Gl) and Z{G*f,_-^) differ is bounded by the right-hand side of (5.7), 
which completes the proof. □ 
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We can now estimate the expectation of Z{Gl.) conditioned on the value 
of XyuiiG'l). We state only the result we need. (See also (5.12). These 
results can be rewritten as estimates of conditional expectations.) 

Lemma 5.4. Suppose ^^d^ = 0{N). With notations as above, 
E{Z{G*)Je,iG*)) 

(I m X 

[N-^ + Y,d-')J{^^u.\{nl^. (5.9) 
veFi 1=1 j=l ^ 

Proof. We can write Xyyj{G*)- = X^ae^ -^o' where A is the set of all ordered 
/c-tuplcs of disjoint pairs (x, y) of half-cdgcs with x belonging to v and y to w, 
and la is the indicator that the k pairs in a all belong to the configuration. 
By symmetry, E(Z(G*) | = l) is the same for all a ^ A; since Gj^ is 
obtained by conditioning G* on I^, for a specific a, we thus have E(Z(G*) | 
7^ = 1) = ^[Z{Gl)) for all a. Consequently, 

E (X,^(G*)^Z(G*)) = E(/«Z(G*)) = E(Z(G*) | = l) P(/. = 1) 

= ^z{Gl)) Y = ^z{Gi)) ^x,^{G*)^). 

(5.10) 

We write the error term on the right-hand side of (5.7) as 0{R). Since 
(5.7) is uniform in k, and Gq = G*, Lemma 5.3 yields 

E{Z{Gl)) = E{Z{G*)) + 0{kR). (5.11) 

We now use Lemma 3.1(ii) and (i) and find, for any j, using (5.10) and 
(5.11), 



E(Z(G*) • l[X,UG*)=j]) = X;(-1)'"'Q)^E(X,^(G*)^Z(G*)) 

k=j 

= f;(-l)'-^' Q y 1E(X„^(G*)^) (E Z{G*) + 0{kR)) 

k=j 

= E{Z{G*))F{X,UG*)=j)+0 QiE(X„^(G*)^)/ci?j . 
By (5.2), the sum inside the last O is at most 
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Since fi^w = 0(1) by (2.3), we thus find, uniformly in j > 1, 

= E(Z(G*)) r{XUG*) = j) + 0{^ii^R/{j - 1)!) , (5.12) 

wliich by summing over j > 2 yields, again using /ly^ = 0(1) and recalling 
that vw = ei, 

K{Z{G*)JeAG*)) = HZiG*))E{JeAG*)) +0{i^lR), 
the sought estimate. □ 

Proof of Lemma 5.1. As said above, we use induction onl + m. The result 
is trivial if / + m = or / + m = 1. If m > 1, we use Lemma 5.4; the 
result follows from (5.9) together with the induction hypothesis applied to 
E(Z(G*)) and the estimate E(Jei(G*)) = 0{^il^) from Lemma 5.2. 

If m = 0, we study a product /^^ of loop indicators only. We then 
modify the proof above, using loops instead of multiple edges in the condi- 
tionings. More precisely we now let G^, be G* conditioned on the configura- 
tion containing the k specific pairs (ti^^*"^), tt^^*)), i = 1, . . . , /c, of half-edges 
at u. We couple Gj^ and as above (with obvious modifications). In this 
case, the switching from Gj^ to Gj^_^ cannot create any new loops. Hence, if 

Z := rTi=2-^«i' have Z(G^) > Z(G^_-|^). We obtain (5.8) exactly as be- 
fore, and since Lemma 5.2 still holds, this shows that Lemma 5.3 holds, now 
with Fi = and the error term 0{N~^ Y[i=2 I* follows that Lemma 5.4 
holds too (with Jg^ replaced by I^i and Fi = 0) by the same proof as above. 
This enables us to complete the induction step in the case m = too. □ 

Remark 5.5. Similar arguments show that Lemmas 5.3 and 5.4, with ob- 
vious modifications, hold in this setting, where we condition on loops at u, 
also for m > 0. A variation of our proof of Lemma 5.1 would be to use 
this as long as / > 0; the result in our Lemma 5.3 then is needed only when 
I = 0, which eliminates cases (ii), (v), (vi), (vii) from the proof. On the 
other hand, we have to consider new cases for the loop version, so the total 
amount of work is about the same. 

Remark 5.6. When conditioning on loops, it is possible to argue directly 
with conditioning on = k, using a coupling similar to the one for 
above; we thus do not need the detour with G^ and Lemma 3.1 used above. 
However, as said above, in order to treat multiple edges, the method used 
here seems to be much simpler. A possible alternative would be to use the 
methods in McKay [10] and McKay and Wormald [11, 13]; we can interpret 
the arguments there as showing that suitable switchings yield an approxi- 
mate, but not exact, coupling when we condition on exact numbers of edges 
in different positions. 

Remark 5.7. A small example that illustrates some of the complications 
when conditioning on a given number of edges between two vertices is ob- 
tained by taking three vertices 1, 2, 3 of degree 2 each. Note that if X12 = 1, 



16 



SVANTE JANSON 



then the multigraph must be a cycle; in particular, = 0. On the other 
hand, = 1 is possible for X12 = 0; this shows that it is impossible to 
couple the multigraphs conditioned on X12 = and on X12 = 1 by moving 
only two edges as in the proof above. Note also that X3 = is possible also 
when X12 = 2; there is thus a surprising non-convexity. 



6. The proofs are completed 

Proof of Theorem 1.4- We begin by observing that the two expressions given 
in (1.4) and (1.5) are equivalent. Indeed, if we define A by (1-6), then 

i i<j i i,j i 

_l^di{di-l) 1 di{di - l)dj{dj - 1) ^ dlidi - 1)2 N 



2 ^ 2iV 4 IP ^ 4Ar2 

^ + ^2 E^d^id^-lr 
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and 



AHA=(a + 1)^-1 = -1. (6.1) 

V 2/44V2Ary4 ^ ^ 

We note for future reference that 

i i i 

and A — log(l + A) > when A > 0, and thus the right-hand side of (1.4) 
can be estimated from above by 

/I Vd^X 

right-hand side of (1.4) < exp - - ^f-^ + o(l). (6.2) 



AN 

Similarly, log(l -|- A) - A -|- ^A^ > when A > 0, and thus 

right-hand side of (1.5) > exp(^-(^^)^^ + o(l). (6.3) 

In particular, since we just have shown that these two right-hand sides are 
the same, they tend to if and only if df /N 00. 

Next, suppose first that J2i di = 0{N). RecaU Y{G*) = E«ey„ + 
J2eeE„ 'Je{G*) defined in (2.1). As said above, Lemma 5.1 shows that the 
random variables luiG*) and Je{G*) are almost independent. We can com- 
pare them with truly independent variables as follows. 

Let lu and Jg be independent 0-1 valued random variables such that 
P(J„ = 1) = P(/„(G*) = 1) and P(Je = 1) = P(Je(G*) = l), and let 

uev„ eeEn 
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Fix A; > 1. We use Lemma 5.1 for all pairs {l,m) with I + m = k and 
sum (5.4) over all such (Z, m) and all distinct ui, . . . ,ui and ei, . . . , em, mul- 
tiplying by the symmetry factor (^) , and noting that the first term on the 
right-hand side of (5.4) can be written Edlj-^wi Tlj <^j)- This gives 

E(y(G*)^) =E(F^) + o(^((iV-i + ^d;^)nMn,n<.)), (6.4) 

^ ^ veF 1=1 j=i ^ ^ 

summing over all such Z,m, (iti)i, (ej)j and with F depending on them as in 
Lemma 5.1. 

Consider one term in the sum in (6.4), write as usual = VjWj, and let 
H be the multigraph with vertices V{H) = {ui} U {vj,Wj} and one loop at 
each Ui and two parallel edges between Vj and wj for each j < m. Let dy-u 
be the degree of vertex v in H and note that each degree d^j-H is even, and 
thus at least 2, and that F is the set of vertices with dy-n > 4. We have 

I m 

ll^'u.ll^^l=^~'^''^ n (6-5) 

i=l j=l veV{H) 

where e{H) = I + 2m is the number of edges in H. 

We group the terms in the sum in (6.4) according to the isomorphism type 
of H. Fix one such type H, and let it have h vertices with degrees ai, . . . , a/j 
(in some order) and b edges; thus b = ^ Si=i % • The corresponding H are 
obtained by selecting vertices vi, . . . ,Vh & Vn, these have to be distinct and 
it may happen that some permutations give the same H, but we ignore this, 
thus overcounting, and obtain from (6.5) that 

Hen \=i j=i ^ v^,...,VHeVr, ^ j=i ^ j=i \jeVn ^ 

= O (^AT-^+S.- = 0(1) (6.6) 

by (2.4), since each aj > 2. 

Furthermore, let G := {i e {1, . . . ,h} : at > 4}. Thus, if H is obtained 
by choosing vertices vi, . . . ,Vh & Vn, then F = {vi : i e G}. Consequently, 

jjc/i vGF 1=1 j=l ^ vi,...,v,^£V„ ieG ^ j=i ^ 

by (2.4), since each > 4 if i € G and thus aj — 5ij > 2 for every j. 

Combining this with (6.6), we see that the sum in (6.4), summing only 
over H E H, is 0{N~^/'^). There is only a finite number of isomorphism 
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types Ti for a given A;, and thus we obtain the same estimate for the full 
sum. Consequently, (6.4) yields 

E(F(G*)^) = E(F-) + 0{N-^/^), (6.7) 

for every fixed k. 

We use Lemma 3.2 with Y and Y{G*') (in this order). We have just 
verified (3.3). To verify (3.2) we take R = 2 (any R < oo would do by a 
similar argument) and find, using (5.3) and (2.2) 



(2^) = n ^2^" n ^^'^^ = n n i^+^^e) 

ueVn e€En ueVn e€E„ 

< Y[ exp(E7„) Y[ exp(E Je) = exp(^ ^ E7„ + ^ E Je j 

s-p(e|+ E #)=o(i). 



Consequently, Lemma 3.2 applies and yields 
F{G* is simple) = F{Y{G*) = 0) = P(F = 0) + o(l) 

= H P(7„ = 0) J] P( Je = 0) + 0(1) 

= expT ^ logP(7„ = 0) + ^ logP(Je = 0)] + 0(1). 
Furthermore, Lemma 4.1 yields 

-[Y, logP(7„ = 0) + ^ logP(Je = 0)") 

= E "^"^IV + E (-log(l + A..) + A„^) 



4Ar 



where the two error terms ore 0{N-^/^) by (2.4). 

This verifies (1.4) and thus Theorem 1.4 in the case '^^df = 0{N). 

Next, suppose that Yli —>■ 00. Fix a number A> 2. For all large N (or 
u in the formulation of Theorem 1.1), > AN, so we may assume this 

inequality. 

Let j be an index with dj > 1. We then may modify the sequence {di)i 
by decreasing dj to dj — 1 and adding a new element dn+i = 1. This means 
that we split one of the vertices in G* into two. Note that this splitting 
increases the number n of vertices, but preserves the number N of edges. 
We can repeat and continue splitting vertices (in arbitrary order) until all 
degrees di < 1; then J2i df = J2idi = 2iV. 
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Let us stop this splitting the first time d? < AN and denote the 
resulting sequence by (di)". Thus < AN. Since we have assumed 

> AN, we have performed at least one split. If the last split was at 
j, the sequence preceding (o?j)" is (dj + ^ij)""^, and thus 

n— 1 h h 

AN < ^{di + 5ijf = ^di + 2dj + 1-1 <J2di + 2VIiV, 

1=1 i=l i=l 

because dj < c/| < AN. Consequently, 

^A^ - 2VAN <^dl < AN 

i 

and thus, in the limit, Yli^^/^ ~^ 

Let G* = G*{h,{di)f). Since = 0{N), we can by the already 

proven part apply (1.4) to G* and obtain, using (6.2), 

P(G* is simple) < exp(^^ - ^) + o(l) = expQ - |) + o(l). 

Furthermore, since G* is constructed from G* by splitting vertices, G* is 
simple whenever G* is, and thus F{G* is simple) < P(G* is simple). Conse- 
quently, 

limsupP(G* is simple) < limsupP(G* is simple) < exp(— (yl — 2)/4). 

Since A can be chosen arbitrarily large, this shows that if Ylii^l/^ ~^ 
then P(G* is simple) 0. 

Combined with (6.2), this shows that (L4) holds when X^jd^/iV — oo 
(with both sides tending to 0). 

Finally, for an arbitrary sequence of sequences (cij)", we can for every sub- 
sequence find a subsequence where either J^i^f/N = 0{1) or Yli^f/N — *■ 
oo, and thus (1.4) holds along the subsequence by one of the two cases treated 
above. It follows by the subsubsequence principle that (1.4) holds. □ 

Proof of Theorem 1.1. Part (ii) follows immediately from Theorem 1.4 and 
(6.2), (6.3). 

If we apply this to subsequences, we see that liminf P(G* is simple) > if 

and only if there is no subsequence along which Ylii^t'^)'^ /^i' ~^ which 
proves (i). □ 

Proof of Corollary 1.5. The two expressions arc equivalent by (6.1). 

If ^idf/N oo, then the right-hand side tends to 0, and so does the 
left-hand side by Theorem 1.1. Hence, by the subsubsequence principle, it 
remains only to show the result assuming (if = 0{N). In this case we 
have 

j:,d^{d, - 1)2 ^ (max,d,)2E.<^ ^ ^^1) 



16A^2 - N N 
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and, since log(l + x) — x + ^x^ = 0(x'^) for a; > 0, 

J](log(l + A,,) - A,,- + iAf,.) = O (E 4) 

i<j ^ i,j 

and 

Hence, in this case, the formula in Corollary 1.5 follows from (1.5). □ 



7. POISSON APPROXIMATION 

As remarked in the introduction, when maxj di = o{N^/^), it is easy to 
prove that (1.2) implies (1.1) by the Poisson approximation method of Bol- 
lobas [2, 3]. Since this is related to the method above, but much simpler, 
and we find it interesting to compare the two methods, we describe this 
method here, thus obtaining an alternative proof of Corollary 1.5. We as- 
sume X^jd? = 0{N) throughout this section. 

The main idea is to study the random variable 

?^:=E^"+E (t)' (7-1) 

which counts the number of loops and pairs of parallel edges (excluding 
loops) in G* (we omit the argument G* in this section). Compare this with 
Y defined in (2.1), and note that 

G* is simple ^ y = ^ y = 0. 

Theorem 7.1. Assume that N 00, Y^^df = 0{N) andm.ax.idi = o(N^/'^). 
Let A := S"=i (2) (1.6). Then 

dTv(>^,Po(A + A2)) ^0, (7.2) 

and thus 

P(G*(n, (di)^) is simple) = ¥{Y = 0) = exp(-A - A^) + o(l). (7.3) 

If A ^ A for some A e [0, 00), then (7.2) is equivalent to Y Po(A+A^). 
(By the subsubsequence principle, it suffices to consider this case.) 

Sketch of proof . We can write Y = J2aeA-^°' X]/3eB'^/3' 'where A is the 
set of all pairs {u^^\u^^^ of half-edges (correponding to loops), and B is 
the set of all pairs of pairs {{v^^\w^^^},{v^''\w^''^}} of distinct half-edges 
(corresponding to pairs of parallel edges). 
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Thus, similarly to (4.1), 

aeA l3eB 

E dujdu - 1) , 1 1 dv{dy - l)dw{dy, - 1) 
^^^2{2N-1) 2^2 (2Ar- l)(2A^-3) 

= A + A2 + o(1). (7.4) 

Moreover, it is easy to compute the expectation of a product of the 
form E(/ai • • • laiJ/Si • • • <^/3m); it is just the probability that a random con- 
figuration contains all pairs occurring in ai, . . . ,ai, (3i, . . . , (3m ■ If two of 
these pairs intersect in exactly one half-edge, the probability is 0; other- 
wise it is (2iV)""^(l -I- 0{1/N)), where b is the number of different pairs. 
(Note that we may have, for example, (3i = {{v^^\w^^^},{v^'^\w^^^}} and 
^2 = {{v^^^ , uM^}, {v^^\ u/'^)}}, with one pair in common; thus b < I + 2m, 
but strict inequality is possible.) 

We can compute factorial moments E,(Y-) by summing such expectations 
of products with I + m = k. For each term E(Iq,^ • • • Jfs^), let H be the 
multigraph, with vertex set a subset of Vn, obtained by joining each pair 
occurring in qi . . . , (taking repeated pairs just once) to an edge, and then 
deleting all unused (i.e., isolated) vertices in Vn- It is easy to estimate the 
sum of these terms for a given H, and we obtain 0(^N~^^^^ Tlv&V{H) d^^'") ^ 
in (6.5). As in the proof of Theorem 1.4, wc then group the terms according 
to the isomorphism type H of H. (There are more types H now, but that 
does not matter.) 

Since now maxj c/j = o{N^/^), (2.4) is improved to 

J24 = oiN'^/') (7.5) 

V 

for every fixed k > 3, and it follows that the sum for a given TC is o(l) as 
soon as Ti. has at least one vertex with degree > 3. The only remaining case 
is when H, and thus H, consists of I and m vertex-disjoint loops and double 
edges; in this case 

I m 

n -^/^O = ((2^ - 1) . . . (2Ar - 2/ - 4m + 1))"' 
i=i j=i 

I m 

= (1 + o(i/iv)) n E n ^ -^Pi ■ c^-e) 

Similarly, we can expand (EY)^ = (Z^Qg^E/o, -|- X]/3gB ^ •-'z?) as a sum 
of terms rTi=i E lat YYjLi ^ with l+m = k. (Now, repetitions are allowed 
among Ui and Pj.) If we introduce H and 7i as above, we see again that 
the sum of all terms with a given H is o(l) except when H consists of I and 
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m vertex-disjoint loops and double edges. The terms occurring in this case 
are the same as in (7.6), and hence their sums differ by 0{1/N) only (since 
these sums are 0(1), see (7.4)). 

Consequently, summing over all 7i and using (7.4), for every A: > 1, 

E(y^) = {V.Yf + o{l) = {K + K^f + o{l). 

If A A, this shows Y — ^ Po(A + A^) by the method of moments. In 
general, we obtain (7.3) and (7.2) by Lemma 3.2 and Remark 3.3. □ 

Remark 7.2. This argument further shows that, asymptotically, the num- 
ber of loops is Po(A) and the number of pairs of double edges is Po(A^), 
with these numbers asymptotically independent. 

In order to compare this method with the one in the preceding sections, 
note that Y >Y and that Y = Y \i and only if there are no double loops or 
triple edges. It is easy to see that if X^jd? = 0{N) and maxjdj = o{N^^^), 
then, using (5.2) and (7.5), 

P(F 7^Y)<Y^ E(Xf ) + ^ E(X|) 

ueVn e&En 

= o(%^)+o(^%Jl!)=o(l), (7.7) 

so in this case the two variables are equivalent asymptotically. In particular. 
Theorem 7.1 is valid for Y too. It is evident that the argument to estimate 
factorial moments of Y in the proof of Theorem 7. 1 is much shorter that the 
argument to estimate factorial moments of Y in the preceding sections. The 
reason for the difference is the ease with which we can compute ^(lai • • • Jpm) 
for a random configuration. Hence the proof of Theorem 7.1 is preferable in 
this case. 

On the other hand, if maxdj = 6(A''^/^), still assuming Y^^d^ = 0{N), 
there are several complications. Let us for simplicity assume that di > 

d2 > . . . , and that di ~ ciN^^ with ci > 0. Then Xi Po(cf /4), so 
limP(X„ > 1) > and (7.7) fails. 
Moreover, cf. (7.4), 

i i<j 

SO we can write (1.4) and (1.5) as 
P(G* is simple) = exp(- + ^(log(l + Xij) - Xij + iAf^)) + o(l). 

Except in the case d2 = o(iV^/^), we cannot ignore the terms with X^j in the 
exponent; if, say, d2 ~ C2N^^^ with C2 > 0, then A12 — > C1C2/2 > 0. Conse- 
quently, Theorem 1.4 shows that in this case, ¥{Y = 0) = P(G* is simple) 
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is not well approximated by exp(— EY), which shows that Y is not asymp- 
totically Poisson distributed. (The reason is terms like ("^2^) C^--'-)' where 

^12 Po(ciC2/2).) 

Further, we have shown in Section 6 that Y asymptotically can be re- 
garded as the sum Y of independent indicators, but in this case limP(/i = 
1) > 0, and thus these indicators do not all have small expectations; hence 
Y is not asymptotically Poisson distributed. 

Any attempt to show Poisson convergence of either y or y is thus doomed 
to fail unless maxj dj = o(iV^/^). It seems difficult to find the asymptotic 
distribution of Y directly; even if we could show that the moments converge, 
the moments grow too rapidly for the method of momemts to be applicable 
(at least with the Carleman criterion, see Section 4.10 in [6]). This is the 
reason for studying Y above; as we have seen above, the distribution is 
asymptotically nice, even if our proof is rather complicated. 

Acknowledgement. I thank Bela Bollobas and Nick Wormald for helpful 
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